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Abstract 

The final oscillation analysis of the complete set of data collected by CHORUS in the 
years 1994-1997 is presented. Reconstruction algorithms of data extracted by electronic 
detectors were improved and the data recorded in the emulsion target were analysed by new 
automated scanning systems, allowing the use of a new method for event reconstruction in 
emulsion. CHORUS has applied these new techniques to the sample of 1996-1997 events 
for which no muons were observed in the electronic detectors. Combining the new sample 
with the data analysed in previous papers, the overall sensitivity of the experiment to the 
v T appearance is thus improved. In a two-neutrino mixing scheme, a 90% C.L. upper limit 
of sin 2 29 ^ T < 4.4 x 10~ 4 is set for large Am 2 , improving by a factor 1.5 the previously 
published CHORUS result. 
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1 Introduction 

The CHORUS experiment was designed to search for — ► v T oscillations through 
the observation of charged-current interactions v T N — > r~X followed by the de- 
cay of the r lepton, directly observed in a nuclear emulsion target. The experiment 
aimed at achieving maximum sensitivity on the effective mixing angle for values of 
the mass parameter Am 2 larger than 10 eV 2 /c 4 . This particular choice was based on 
the hypothesis that the neutrino mass could contribute to the solution of the Dark 
Matter puzzle H1I2II . A short-baseline experiment in the CERN SPS Wide Band 
Neutrino Beam [O was well suited for this search. Another experiment, NOMAD, 
used the same beam and searched for v T appearance using a purely electronic tech- 
nique flU. 

The CHORUS experiment took data from 1994 to 1997. A first phase of data anal- 
ysis ('Phase F) was performed and no evidence for oscillations was found H5I6I7I . 
Owing to several improvements in automated emulsion scanning and in the event 
recontruction, it was considered worthwhile to perform a new and more complete 
analysis ('Phase II') of the data collected in the 1996-1997 period. 

In this paper, we report the analysis performed on the events for which no muons 
were observed in the electronic detector (0/x events). In particular, a dedicated 
search for decays into three charged hadrons was also performed. Combining this 
new sample with the data samples analysed and published in previous papers, i.e. 
with a muon observed in the electronic detector (1//) for the whole period of data 
taking and a small 0/x sample collected in 1994-1995, leads to an improvement of 
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the CHORUS sensitivity to — > z/ r appearance by a factor 1.5 and to u e — > z/ r 
appearance by a factor 1.2. 

Although it is now established that — > i/ T oscillations occur at Am 2 ~ 10~ 3 , this 
paper intends to show the capabilities of an hybrid emulsion experiment and that 
the goal for which CHORUS was designed has been reached. 



2 The experimental setup 

The CERN SPS Wide Band Neutrino Beam consists essentially of with a con- 
tamination of about 5.1% V^, 0.8% of u e , and 0.2% of V e , while the flux of v T is 
negligible [ 8]. The average neutrino energy is 26 GeV, well above the r production 
threshold. 

The CHORUS detector is a hybrid setup which combines a nuclear emulsion target 
with various electronic detectors such as trigger hodoscopes, a scintillating fibre 
tracker system, a hadron spectrometer, electromagnetic and hadronic calorimeters, 
and a muon spectrometer 

Thanks to its micrometric granularity, nuclear emulsion allows a precise three- 
dimensional reconstruction of the neutrino interaction vertex as well as of the decay 
vertices of associated short-lived particles. It is a powerful detector technique to di- 
rectly observe r decays which, in this experiment, occur on average at a distance 
of 1.7 mm from the interaction vertex. 

The scintillating fibre tracker system consists of two components: the target tracker 
(TT), which is the major tool for locating in the emulsion the region where a neu- 
trino interaction has occurred, and the diamond tracker (DT), a set of trackers as- 
sociated with the hadron spectrometer. 

The hadron spectrometer measures the charge and momentum of charged parti- 
cles using an air-core magnet. The calorimeter is used to determine the energy of 
electromagnetic and hadronic showers. The magnetized iron muon spectrometer 
determines the charge and momentum of muons. 

The emulsion target has an overall mass of 770 kg and is segmented into four 
stacks. Each stack consists of eight modules with 36 plates of size 36 cm x 71 cm. 
Each plate has a 90 yum plastic base coated on both sides with a 350 /im emulsion 
layer. Each stack is followed by three interface emulsion sheets having a 100 /xm 
emulsion layer on both sides of a 800 thick plastic base, and by a set of TT 
planes. The TT predicts particle trajectories on the interface emulsion sheets with 
a precision of about 150 fim in position and 2 mrad in track angle. 

The emulsion scanning is performed by computer-controlled, fully automated mi- 
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croscope stages equipped with a CCD camera and a read-out system called 'track 
selector' H10I111I . The last generation of automated microscopes used in CHORUS 
experiment is the UTS ('Ultra Track Selector') [12J. In order to recognize track 
segments in the emulsion, a series of tomographic images is taken by focusing 
at different depths in the emulsion layer. The digitized images at different levels 
are shifted according to the predicted track angle and then added. The presence of 
aligned grains forming a track is detected as a local peak in the grey-level of the 
summed image. The track-finding efficiency is higher than 98% for track slopes 
less than 400 mrad IfTTTl. 



3 The event reconstruction 



The event reconstruction algorithms have been improved significantly with respect 
to the Phase I analysis described in detail in Ref. 0. Here, we only summarize the 
essential points. 

The event reconstruction starts with the pattern recognition in the electronic detec- 
tors. Tracks are found in the TT, in the calorimeter and in the muon spectrometer. If 
a muon is detected in the downstream detectors the event is classified as other- 
wise it is called O/i. Vertices are reconstructed by the electronic detectors using the 
points of closest approach of the tracks in the TT immediately downstream from the 
target. The main vertex is the most upstream one and is selected if it contains either 



a muon or a hadron surviving criteria which differed in the two phases UU . Such 
particles or the muon are used as so-called 'scan-back' tracks. The impact points 
of the scan-back tracks are predicted on the most downstream interface emulsion 
sheets, to initiate their follow back in emulsion. 

As in previous analyses the event reconstruction in emulsion starts with a procedure 
called 'vertex location' . As a first step, all the tracks in the interface emulsion sheets 
within an area of 1 mm 2 centred on each scan-back track prediction are collected. 

The correct association of emulsion tracks with TT tracks, both in position and 
direction, requires a precise alignment of the interface emulsions with respect to 
the fibre-trackers. The alignment parameters are obtained by finding the best match 
between the full set of TT predicted tracks and the full set of collected emulsion 
tracks. Once an emulsion track is associated to a scan-back track, mainly on the 
basis of their common direction, it is followed upstream from one emulsion plate 
to the next within a greatly reduced scanning area as the resolution improves. In the 
emulsion stacks, the scanning area reduces to a square of 50 fim side. The plate- 



14 In Phase I, Ofi events, at least one hadron with a momentum less than 20 GeV/c was 
required. In Phase II the most isolated hadron was selected regardless of its charge and 
momentum. 
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to-plate alignment is first obtained by a coarse adjustment of reference marks, and, 
as the scanning procedure continues, is refined by aligning track maps measured 
in subsequent plates. For this alignment we used tracks of particles generated by 
neutrino interactions during the whole data taking, of cosmic rays and of muons 
coming from neighbouring beams. The scan-back procedure stops when a searched 
track is not found in two consecutive plates. The most downstream one is defined 
as the "vertex plate". This procedure for locating the vertex has an efficiency of 
34% (Table [[])• 

Once the vertex plate is found, a new scanning technique is applied on the 1996— 
1997 data sample. Originally developed for the DONUT experiment |[T6ll , this tech- 
nique, called 'NetScan' IfPTll , is described in detail in Ref. H13I18H together with the 
event location procedure. Its application to the 0/x sample is also described in the 
following section Uf. . 

3. 1 The NetScan technique 

Once the vertex plate is identified, the UTS performs a scan of the emulsion vol- 
ume around the vertex position, recording, for each event, all track segments within 
400 mrad with respect to the orthogonal direction of the plates. In each plate only 
the most upstream 100 impart is scanned. The scanned volume is 1.5 mmx 1.5 mm 
wide and 6.3 mm long in the beam direction, corresponding to eight emulsion 
plates. This volume contains the vertex plate itself, the plate immediately upstream, 
and the six plates downstream the vertex plate. The plate upstream the vertex acts 
as a veto for passing through tracks. The six plates downstream from the vertex act 
as decay volume and are used to detect the tracks of the decay daughters. The scan 
area is centred on the scan-back track stopping point. 

The number of track segments (coming from particles induced by neutrino inter- 
actions, cosmic ray particles and muons from neighbouring beams) found in each 
plate depends on the position of the NetScan volume with respect to the beam cen- 
tre and on average is 920. 

The task of the NetScan event reconstruction is to select the segments belonging to 
the neutrino interaction under study, out of this large number of background track 
segments. To this end, a first plate-to-plate alignment is performed by compar- 
ing the pattern of segments in a plate with the corresponding pattern in the next 
upstream plate. Each segment found in one plate is extrapolated to the next plate 
where a matching segment is looked for within about 4 /xm in position (3a of align- 
ment resolution) and 20 mrad in angle. If none is found, the straight-line extrapo- 
lation is tried one plate further upstream. 

15 The NetScan technique has been used also for the lfi data sample, mainly to search for 
charmed particles decays 1141151 . 
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A second and more accurate inter-plate alignment is performed using tracks passing 
through the entire volume after the connection of all matched segments. These 
tracks come from muons associated with the neutrino beam or with charged particle 
beams in the same experimental area. After this fine alignment, the distribution of 
the residuals of the segment positions with respect to the fitted track has a standard 
deviation width of about 0.45 /mi. 

After rejection of isolated track segments, typically about 400 tracks remain in the 
volume. The majority of these are tracks (mainly Compton electrons and 5-rays) 
with a momentum less than 100 MeV/c. These background tracks are rejected on 
the basis of the x 2 of a straight-line fit to the track segments. The final step is the 
rejection of the tracks not originating from the scan volume. After this filtering, the 
average number of remaining tracks is about 40. The number of connected tracks 
and the position and slope residuals provide a good measurement of the alignment 
quality. About 84% of the events passes the quality cuts ('NetScan acquisition ac- 
cepted'). 

The reconstruction algorithm then tries to associate the tracks to vertices with a 
'pair-based' method: after selecting pairs of tracks having minimal distance less 
than 10 /im, a clustering among them is performed and vertex points are defined. 
After this clustering, a track is attached to a vertex if its distance from the vertex 
point (called hereafter 'impact parameter') is less than 3 /im. At the end of the 
procedure, one defines a primary vertex (and its associated tracks) and in case one 
or more secondary vertices to which 'daughter tracks' are attached. 

About 77% of the events, where the NetScan acquisition was accepted, have a ver- 
tex reconstructed in the NetScan volume. The loss is accounted for by the fact 
that a background track might have been selected for the scan-back procedure. The 
numbers of events at the various stages of the procedure are given in Tabled! 

Table 1 

Results of the reconstruction of the Ofi sample 

Stage of reconstruction number of events 

Interface emulsion scanned 102544 

Vertex plate found 35039 

NetScan acquisition accepted 29404 

Vertex reconstructed 22661 

3.2 The decay search 

Once tracks and vertices are reconstructed in the NetScan volume, a search for 
decay topologies is performed. In Fig.[T]a sketch of a v T Monte Carlo event inside 
the NetScan volume is shown: a v T (not drawn) interacts producing a r lepton 
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and other particles. After three plates, the r lepton decays producing one charged 
particle (kink topology, i.e. an observed abrupt change of direction in the track). 




1500flm 



8 7; vertex 6 5 4 3 2 1 1500m. m 



lOOum scanned 




Fig. 1. The NetScan volume and a Monte Carlo event with one-prong (i.e. kink) r decay. 

In the decay search an emulsion track is called a 'TT-confirmed track' if its trajec- 
tory matches that of a TT reconstructed track, or if it is connected to a sequence of 
TT hits. A reconstructed vertex is a 'TT-confirmed vertex' if it contains at least one 
TT-confirmed track. 

The decay search is done in two stages: 

i) an automated search on the whole sample; 

ii) a computer-assisted visual inspection on the sample selected in the previous 
stage. 

After the automated search, 754 events are selected as belonging to one of three 
categories: 

• Short decays. Events for which there is only one TT-confirmed vertex in the 
NetScan volume and at least one isolated TT-confirmed track having an impact 
parameter with the primary vertex inside the range [7.2, 100] fim. 

• Long decays. Events with a one-prong decay topology ('kink') directly observed 
in the NetScan volume. This occurs when the parent track r is observed in at least 
one emulsion plate. The daughter track must be TT-confirmed and its minimal 
distance to the parent track is required smaller than 6 fim. In case a primary 
vertex is reconstructed, it is required to be TT-confirmed and the parent track 
must have an impact parameter less than 3 /im with this vertex. 

• Multiprong decays. Events in which NetScan reconstructs at least two vertices, 
the secondary vertices being candidates for multiprong r decays. The quality of 
the vertices is ranked according to the number of TT-confirmed tracks originating 
from them and it is used as a parameter for the selection of multiprong decay 
topology. The parent angle, in the case of a neutral particle decay deduced from 
the line connecting the primary and secondary vertices, should be within 400 
mrad from the beam direction. 
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The computer-assisted visual inspection aims at clearly establishing the topology 
of the primary and secondary vertices. To avoid ambiguous topologies, the flight 
length of the parent track is required to be larger than 25 /im. For vertices which 
are not both TT-confirmed, a more stringent cut, which is a function of the flight 
length of the assumed parent, is applied to remove the random association of two 
unrelated vertices. 

An event is selected if one of the following topologies is detected: 

• CI : a track with a kink (interaction or decay) of at least 50 mrad 

• C3 : a track with a three-prong (interaction or decay) topology 

• V2 : a two-prong neutral interaction (or decay) is detected 

The numbers of events observed in these topologies are shown in Tabled Decays 
of the r lepton are characterized by either one or three charged particles in the 
final state, CI or C3. The events with a V2 topology will be used for background 
calculations. 

Table 2 

0// events selected after visual inspection. 

Topology number of events 

CI 59 
C3 48 
V2 99 

3.3 Background evaluation 

Physical processes that can mimic a r decay are: 

• decays of charmed charged particles like D + , and A+, if the primary pT is 
not identified and the charge of the charmed particle is not determined; 

• white interactions, i.e. interactions of hadrons without any other visible activity 
at the vertex (recoil or Auger electrons); 

• E decays which affect only the CI topology. 

The charmed particle contribution to the background has been evaluated from pro- 
duction cross-section of D + , and A+ in neutrino interactions and their decay 
branching fractions into one or three charged particles |fi~9)l . The rate of D° produc- 
tion in neutrino charged-current interactions has been measured in CHORUS [14J. 
Using as the normalization the observed number of events with a V2 topology given 
in Table [2l the number of one and three-prong decays from charm particles in the 
O/i sample is given by 
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a D o ev2 

where i = 1 or 3, (Jq/o"do = 1.03 ± 0.16 is the ratio of charged charm and D° 
cross-sections lfi~9l and e V 2, eci and e C 3 are the efficiencies in detecting the specific 
topology, including the branching ratios, computed by a Monte Carlo simulation: 

eci/ev2 = 0.28 ± 0.08, e C3 /e V 2 = 0.45 ± 0.09. 

The number of V2 events has been corrected to account for K° and A decays lfi~8l . 
This yields a prediction of JV C1 = 27.2 ± 8.6 and Acs = 44 ± 11. 

One-prong interactions without visible activity ('white kinks', WK) were gener- 
ated and processed through the full simulation chain assuming a hadron interaction 
length A = 24 m EOl . In the lji sample, the ratio of three-prong to one-prong in- 
teractions, both with visible activity, was measured to be ~ 11%. This fraction was 
assumed to be the same also in absence of recoil or Auger electrons [Tl8l . 



The results are summarized in Table |3] 
Table 3 

Expected background and observed data for CI and C3 topologies in the Ofi sample. 



Topology 



0/iCl 
0//C3 



Charm White interactions Other decays Total BG Data 



27.2 ±8.6 24.9 ±2.7 1.11 ±0.34 53.2 ±9.0 59 
44 ±11 2.7 ±0.3 - 47 ±11 48 



3.4 Maximum number of detectable r 



We define as the number of r events that would be observed if their detection 

efficiency and the oscillation probability would be equal to 1. It can be evaluated 
as: 



tfMAX = f $ {Ev) . {Ev) dE = ^.pl 

J ecc 

where $ y (E v ) is the neutrino flux, N cc (93807) is the number of charged-current 
(CC) neutrino interaction events (used only for the flux determination) detected 
in the 1996-1997 run for which NetScan acquisition is completed, and e cc their 
detection efficiency. (c 1/T )/(cr 1/ ) is the ratio between the average total and v T 
CC cross-sections. Efficiencies have been calculated by Monte Carlo. We define 
e loc as the probability to locate a vertex and e sel as the reconstruction and selection 
efficiencies for the various decay channels. For each topology, the number j\T™ ax of 
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r events that would be observed if the oscillation probability would be equal to 1 
is: 

a rmax atMAX \ " tj D 4oc ,_sel 

1\ T = 1\ T ■ oti m ■ e m ■ e m , 

m 

where the index m runs over the various r decay modes contributing to that topol- 
ogy. Branching ratios, efficencies and jV™ ax are given in Table HI Note that muonic 
decays with a misidentified muon also contribute to the CI topology. 

Table 4 

r branching ratios, efficiencies and maximum detectable r events for CI and C3 topologies 
in the 0/x sample. 



Topology 


Channel 


5i? (%) 




jvrmax 

T 




T - 




49.5 


32.6 ±0.5 23.7 ±0.6 


7227 




0/iCl 


T - 


-> e 


17.8 


26.9 ±0.6 21.2 ±1.3 


1922 


9621 




T - 


-> fi 


17.4 


5.9 ±0.3 24 ±3 


472 




0/xC3 


T - 


-> 3h 


15.2 


35.8 ±0.8 43.1 ± 1.5 


4443 


4443 



3.5 Post-scan selection and final data sample 

The signal to background ratio, hence the sensitivity to oscillation, can be improved 
if, in addition to the topology, other information is used, such as: 

• the angle between the parent track and the estimated direction of the hadronic 
shower 

• the length of the parent track 

• the mean angle of the daughter tracks with respect to the parent track 

• when available, the momentum and charge of a daughter particle 

The purpose of the 'post-scan' analysis is to maximize the sensitivity to oscillation 
by an optimum combination of this information. The optimization was performed 
on Monte Carlo simulated samples of signal and background events |fi~8l and led to 
the subdivision of the samples into different subclasses, depending on the observed 
topology. 

3.5.1 The CI topology 

For the CI topology, the transverse momentum with respect to the parent direction 
is useful for the rejection of WK background. 
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Two methods are available for momentum measurements: either with the DT of 
the hadron spectrometer or from the multiple Coulomb scattering (MCS) in the 
emulsion. The DT also measures the particle charge. The MCS method is used if 
no DT measurement is available and if the number of emulsion plates for the MCS 
measurement is greater than 6. The CI sample has thus been divided into three 
subclasses: 

• momentum measurement from DT (P DT ) which also yields a charge measure- 
ment of the sign of the charge which is important for the charm background 
rejection 

• momentum measurement from MCS (P MCS ) 

• no momentum measurement 

The first category has the highest sensitivity to r search and the strongest back- 
ground rejection. The events in which the daughter track has a positive charge are 
rejected. When the momentum is determined, a cut is applied on the transverse mo- 
mentum (Pt > 250 MeV/c) of the daughter kink. A momentum-dependent cut on 
the flight length (Lp) is also applied. Fig.[2]shows a comparison between signal and 
background in the selected (I) and rejected (II) regions in the L F versus P plane. 
All decay vertices beyond 3 mm or with a parent momentum below 2 GeV/c were 
rejected because white kink interactions are likely to populate these regions. For a 
momentum between 2 GeV/c and 5 GeV/c step functions were used for simplicity. 
The effects of this selection on the background rejection and r efficiency are shown 
in Table [5J 



Table 5 

Effect of the selection in flight length-momentum plane on the signal and background, 
without considering any other cut. 



Physical process 


After rejection of events 


WK 


13% 


charm 


34% 


r 


52% 



In addition, the samples have been further subdivided into three ranges, where 
A(p is defined as the angle in the transverse plane between the parent track and 
the mean direction of the other primary tracks after exclusion of the one which 
makes the largest angle with the parent track (the track with a cross in Fig. [3]). A 
comparison between the distributions of A0 in charm and r samples is shown in 
Fig. SI The r signal is characterized by a larger mean value of A0. 
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3.5.2 The C3 topology 



For the C3 topology, the fraction of events for which the three daughter particles 
have their momentum measured is too small to lead to a useful discrimination. 
Instead, use is made of the fact that the mean angle (9) of the daughter tracks 
relative to the parent is strongly correlated to the inverse of the Lorentz boost factor 
7. Thus, the product of the parent flight length L F and (9) is a good parameter 
to discriminate particles of different lifetimes (or apparent lifetime in the case of 
white interactions). The C3 sample has been divided into "short cr" and "long cr" 
subsamples, separated by a cut value L-p ■ (9) = 75 fim. Fig. \5\ shows a comparison 
of the charm and r samples: the short cr sample is enriched in r decays. These 
categories are further divided as above according to ranges. The sample most 
sensitive to a r signal corresponds to events with short cr and large A0. 

Once the different categories are defined the same subdivision is applied to the 
candidate events. This is done according to the 'blind analysis' prescriptions, in 
order to minimize the possibility of biases in the experimental result. 

Table [6] shows the number of simulated and observed events in CI and C3 topolo- 
gies of the new 0/z sample, together with data from the Phase I analysis which 
consist of two samples 0: 

• 1/i events from the 1994-1997 data taking. The background reduction is ob- 
tained by applying a cut on the transverse momentum (P T > 250 MeV/c) of 
the kink daughter. The kink must occur within five plates downstream from the 
neutrino interaction vertex plate. The estimated background is 0.1 events and 
the maximum number of v T events is N™ ax = 5014 (as reported in Table 2 of 
[7J). The background is much smaller than for the O/i sample because of the low 
probability of a wrong measurement of the charge in the muon spectrometer. No 
candidate is observed. 

• 0/i one-prong events from the 1994-1995 data taking. Table 2 of gives the 
background and the maximum number of v T oscillated events for the whole Ofi 
1994-1997 sample. The 1996-1997 data has been reanalysed, as reported in this 
paper. The 1994-1995 data correspond to an estimated background of 0.3 events 
and a maximum number of v T oscillated events of N™ ax = 526. No candidate is 
observed. 



The errors on the expected background are evaluated by taking into account the 
limited statistics of the Monte Carlo sample and the errors on cross-sections and 
branching ratios. 
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4 Results 



From the number of observed candidates, the expected background and the number 
of signal events expected for full oscillation it is possible to compute, for each 
subsample of Tabled the 90% C.L. upper limit on the v T appearance probability 
using a frequentist statistical approach, the so-called Feldman and Cousins unified 
approach EH. 

To evaluate the sensitivity of the experiment, we have simulated a large number of 
experiments taking into account Poisson fluctuations of the expected background. 
The sensitivity is then defined as the average of the 90% C.L. upper limits on the 
appearance probability of all these "experiments". 

The last two columns in Table [6] show the sensitivities S^ T of each subsample and 
an index in decreasing order of sensitivity. The subsamples are then combined, 
with the same prescriptions used by the NOMAD Collaboration |@]|. After including 
the 1 1th most sensitive subsample, the global sensitivity changes only at the percent 
level. Similarly for the u e — ► v T search. In evaluating the global sensitivity, we thus 
consider only the 1 1 most sensitive subsamples for both searches. 

4.1 z/ M — > v T oscillation 

Fig. [6] (left) shows the distribution of the upper limits (at 90% C.L.) obtained by 
simulating 800 experiments. The average value, S = 2.4 x 10~ 4 , corresponds to 
the CHORUS sensitivity to v T appearance. 

The 90% C.L. upper limit on the appearance probability obtained from the data is 

P(z/ M -> v T ) < 2.2 x 1(T 4 , 

indicated by a vertical line in Fig. [6] It is in agreement with the expected sensitiv- 
ity. In fact, in absence of signal events the probability to obtain an upper limit of 

2.2 x 10 -4 or lower is 48%, which means that the number of observed events is 
compatible with the estimated background. The above result improves by a factor 
1.5 the previously published limit P < 3.4 x 10 -4 Q. 

In a two-neutrino formalism, the above upper limit corresponds to the exclusion 
region in the (Am 2 , sin 2 29^ T ) oscillation parameter plane shown in Fig|7](left) and 
to the limit sin 2 26^ T < 4.4 x 10~ 4 for large Am 2 . 

The upper limit obtained by NOMAD [4 J is more stringent than the CHORUS one 
at large Am 2 , whereas for Am 2 lower than 70 eV 2 /c 4 the upper limit is improved 
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by CHORUS owing to its higher efficiencies at low neutrino energies. 



4.2 v e — > u T oscillation 

The SPS neutrino beam contains a 0.8% v e contamination. Assuming that all ob- 
servable v T would originate from this contamination, the above result translates into 
a limit on the v e — > v T appearance probability. The difference in energy between 
the ((E v ) ~ 26 GeV) and u e ((E Ve ) ~ 42 GeV) components leads to a different 
shape of the exclusion plot in the oscillation parameter plane. 

Fig. [6] (right) shows the upper limit distribution obtained by simulating 800 exper- 
iments. The average value, S = 2.5 x 1CT 2 , corresponds to the sensitivity to the 
v e — > v T appearance probability. 

The 90% C.L. upper limit on the appearance probability obtained from the data is 

P{v e -> v r ) < 2.2 x 1(T 2 , 

indicated by a vertical line in Fig. [6] In absence of signal events, the probability 
to obtain an upper limit of 2.2 x 1CT 2 or lower is 49%. The above result improves 
by a factor 1.2 the previously published limit [7 ]. The improvement is smaller than 
in the — > v T case for two reasons. The first one is the variation of the expected 
v e contamination provided by a new simulation of the beam (reduced by a factor 
0.85) [|23l . Furthermore, in the — > v T case the efficiency improved predomi- 
nantly at low neutrino energies and the harder v e spectrum makes this improvement 
less effective. This reason also explains the large difference in the limit relative to 
NOMAD. In a two-neutrino formalism, the above upper limit corresponds to the 
exclusion region in the (Am 2 , sin 2 29 CT ) oscillation parameter plane shown in FigJT] 
(right) and to the limit sin 2 29 eT < 4.4 x 10~ 2 for large Am 2 . 



5 Conclusions 



The final result of the search for v T appearance with the CHORUS experiment has 
been presented. A 90% C.L. upper limit of 2.2 x 1CT 4 was set for the appearance 
probability, improving by a factor 1.5 the previously published CHORUS result. 
In a two-neutrino mixing scheme, this result corresponds to the limit sin 2 29^ T < 
4.4 x 10~ 4 for large Am 2 . With respect to previous CHORUS results, we have also 
improved by a factor 1.2 the upper limit on the (y e — > v T ) oscillation probability. 
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Fig. 2. The flight length Lp versus daughter particle momentum computed for charm (left), 
white kink (centre) and r events (right). The momentum is evaluated by DT (top) or MCS 
(bottom). Only events in region I are selected. 
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Table 6 

The final CHORUS data sample. The first two rows refer to the Phase I analysis, namely 
to the 1/i channel of the whole data taking (1994-1997) and to the Ofi sample collected in 
1994-1995. The new sample, consisting of the Ofi data collected in 1996-1997, is divided 
in CI and C3 topologies which are further divided in subsamples, as described in Section 
13.51 For each subsample, the following quantities are shown: the expected background; the 
maximum detectable number of r events, N? T and N!f T respectively from the and v e 
beam components; the number of data events. The numbers in the last two columns are 
evaluated for the — ► v T search (see Section[4]) and give the sensitivity S^ T of each single 
channel (times 10 4 ) and an index which sorts the sensitivities in decreasing order. 



Category A(j)(rad) 


Background 


Nr 


N? 


Data 


SflT 






t -> 1/x [1994 - 1997 data taking] 


0.100 ±0.025 


5014 


55.8 





4.9 


1 


r -> Ofj, CI [1994 - 1995 data taking] 


0.300 ± 0.075 


526 


5.85 





48 


10 



r -> 0/i CI [1996 - 1997 data taking] 


53.2 ±9.0 


9621 


76.9 


59 








[0;tt/2] 


18.0 ±2.3 


769 


7.67 


26 


211 


15 


No momentum measured 


[7r/2;3vr/4] 


5.00 ±0.73 


708 


5.34 


10 


66 


12 




[37r/4;vr] 


6.2 ± 1.4 


1406 


13.9 


7 


38 


6 


Only MCS momentum measured: 


[0;tt/2] 


4.6 ± 1.1 


991 


6.75 


2 


45 


9 


pMCS > 250 MeV/c 


[7r/2;37r/4] 


1.20 ±0.40 


749 


5.70 


2 


34 


5 


and (pMCS vs Lp) cut 


[37r/4;vr] 


3.3 ± 1.0 


1649 


12.8 


3 


26 


4 


DT momentum measured (Charge -): 


[0;tt/2] 


0.383 ± 0.071 


546 


3.62 





41 


8 


P^ T > 250 MeV/c 


[7r/2;37r/4] 


0.087 ±0.033 


556 


4.24 





38 


7 


and {P DT vs. L F ) cut 


[37r/4;vr] 


0.055 ±0.012 


1023 


7.24 





22 


2 


Pt < 250 MeV/c or Charge + or in region II 


14.6 ± 1.6 


1224 


9.77 


9 







r -» 0^ C3 [1996 - 1997 data taking] 


47 ± 11 


4443 


35.5 


48 








[0;tt/2] 


14.8 ±5.0 


792 


6.63 


17 


133 


14 


Short cr (< 75 fj,m) 


[7r/2;3vr/4] 


6.4 ±3.0 


782 


6.12 


6 


89 


13 




[3vr/4;^] 


1.5 ± 1.5 


1554 


13.4 


4 


23 


3 




[0;tt/2] 


15.5 ±5.0 


386 


2.37 


8 


268 


17 


Long cr (> 75 /xm) 


[7r/2;3vr/4] 


9.8 ±3.9 


336 


2.40 


8 


237 


16 




[37r/4;vr] 


1.7 ± 1.5 


593 


4.62 


5 


62 


11 
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Fig. 3. Definition of the angle A<fi. The dashed line is the mean direction of the primary 
tracks, with the exclusion of the one which has the largest <\> angle relative to the parent 
particle. The A0 distribution is different for charm, white interaction or r events. 
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Fig. 4. Acf) distributions for simulated 0/x CI events: comparison between v T (solid line) 
and charmed events (dashed line). The area is normalized to 1. The vertical lines show the 
applied cut dividing the events in three subsamples. 
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Fig. 5. Lp ■ (9) distribution for simulated Ofj, C3 events: comparison between r (solid line) 
and charm (dashed line) decays. The area is normalized to 1. The vertical line shows the 
applied cut Lp ■ (6) = 75 ^m dividing the events in two subsamples. 




Upper limit at 90% C.L. (10" 4 ) Upper limit at 90% C.L. (10" 2 ) 



Fig. 6. Upper limits obtained at 90% C.L. ll2TTl . in absence of signal events, for 800 simu- 
lated experiments with the CHORUS expected background. The average value corresponds 
to the sensitivity to — > v T (left) and v e — > v T (right) appearance probability. The vertical 
line is the 90% C.L. upper limit obtained from the data. 
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Fig. 7. The CHORUS upper limit on — * u T (left) and v e — > v T (right) oscillation repre- 
sented in an exclusion plot in the oscillation parameter plane. CHORUS results are shown 
as solid lines and are compared with the last results of NOMAD [4] and CHOOZ |[22Tl . 
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